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ABSTRACT 

We have used Zeeman-Doppler maps of the surface field of the young, rapid rotator AB 
Dor (P ro t = 0.514 days) to extrapolate the coronal field, assuming it to be potential. 
We find that the topology of the large-scale field is very similar in all three years for 
which we have images. The corona divides cleanly into regions of open and closed field. 
The open field originates in two mid-latitude regions of opposite polarity separated by 
about 180° of longitude. The closed field region forms a torus extending almost over 
each pole, with an axis that runs through these two longitudes. We have investigated 
the effect on the global topology of different forms of flux in the unobservable hemi- 
sphere and in the dark polar spot where the Zeeman signal is suppressed. The flux 
distribution in the unobservable hemisphere affects only the low latitude topology, 
whereas the imposition of a unidirectional polar field forces the polar cap to be open. 
This contradicts observations that suggest that the closed field corona extends to high 
latitudes and leads us to propose that the polar cap may be composed of multipolar 
regions. 

Key words: stars: activity - stars: imaging - stars: individual: AB Dor - stars: 
coronae - stars: spots 



1 INTRODUCTION 

AB Dor is one of the most comprehensively observed of the 
known young rapid rotators. Long term studies of it s photo- 
metric as w ell as its X-ray variability are available ( jAmadc 
et al. 2001 ). It shows a small (5-13 %) rota tional modulation 
in its X-ray emission ( Kiirster et al. 1997) and a large emis- 
sion measure of 10 52_£,3 cm _3 (Vilhu et al. 2001) consistent 



with a very extended or a very dense corona. The obser- 
vation of large coronal prominences trapped in co-rotation 
with the star between 3 and 5 stellar radii from the rotation 
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axis (Collier Cameron fc Robinson 1989a; Collier Cameron 



et al. 1999) suggest that the corona maintains a complex 
structure out to large distances. 

Several observations also suggest the presence of closed 
loops at high l atitudes. Radio observations taken over a 6 
month period ( Lim et al. 199^ ) have shown two prominent 
peaks in the emission, separated by 180° of longitude. Lim 
et al suggest a model for directed radio emission orginating 
from latitudes around 60°. X-ray emission from such high 
latitudes is also believed to be responsible f or the flare ob- 
served with BeppoSAX ( Vlaggio et al. 2000| ) which showed 



no rotational modulation although the observations spanned 
more than a whole rotation period. Since the loop height de- 
rived from modelling the flare decay phase was only 0.3R*, 
Maggio et al claim that this flaring loop structure must have 
been situated above 60° latitude where it would remain in 
view throughout a rotation cycle. 

These indicators of magnetic loops at high latitude are 
consistent with Doppler images of AB Dor (and many other 
stars - see Strassmeier 1996) which show dark spots at or 
near the pole in addition to spots at low latitudes. Indeed, 
Zeeman-Doppler images (see Fig. |l|) show flux at all latitudes 
on AB Dor (Donati fa Collier Cameron 1997, Donati et al. 
1999) except close to the pole where the Zeeman signal is 
suppressed due to the low surface brightness there. What 
the Zeeman-Doppler images do show at the boundary of the 
dark polar cap is a ring of unidirectional azimuthal field. 
Pointer et al. (2002) suggest that this may be the result of 
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the star's differential rotation dragging out meridional field 
lines at the edge of the polar cap to form an azimuthal ring. 
In this case, the polarity of the field in the azimuthal ring 
would depend on the polarity of the field in the polar cap. 
In particular, the fact that this ring is of uniform polarity, 
would suggest that the field in the polar cap is also of one 
polarity. There is one significant problem with this scenario, 
however. If we place a unipolar field region at the pole, then 
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Figure 1. Zeeman-Doppler images of the surface radial field of AB Dor on 1995 Dec 11-13 (top) and 1996 Dec 23-25 (bottom). The 
scale bar on the right is in Gauss. Note that because the star is inclined at 60° to the observer, there is limited information in the lower 
hemisphere. 
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Figure 2. The global topology of the magnetic field of AB Dor for t he 19 95 d ata s et. The top three panels show a view from longitude 
18° while the bottom panels show a view from longitude 288°. Panels f2(a)| and 2(d) show the separator surfaces that separate the closed 



and open field re gions . No te th at the gap in the surface at the pole is due simply to problems in tracin g fiel d lin es th at pass directly over 
the pole. Panels 2(b) and 2(e) show th e clo sed field lines just inside the separator surface and panels [2(c)| and 2(f) show the open field 
lines just outside it. In panels 2(c) and 2(f) the surface radial map has been painted onto the stellar surface. 



some fraction of those field lines will be forced open by the 
pressure of the plasma they contain. If this "polar hole" 
extends down too far in latitude, then it will not be possible 
to explain the BeppoSAX observations that suggest that the 
closed corona extends to latitudes above 60°. 

In this paper we use the radial magnetic field maps 
shown in Fig. hi to extrapolate the coronal field. Our aim is 



to study the topology of the large scale field and to examine 
whether it is possible to reconcile the observations indicating 
a closed corona at high latitudes with the presence of a large 
polar spot. 
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Figure 3. Images of the combined radial field with the 1995 map in the upper hemisphere and the 1996 map in the lower hemisphere. 
The top two panels show the case where the two maps are aligned such that both hemispheres are symmetric about the equator (positive 
polarity open region at the same longitude in each hemisphere). The lower two panels show an antisymmetric alignment (positive polarity 
open regions 180 degrees apart in longitude in the two hemispheres). Left-hand side panels show the surface radial field (greyscale extends 
from black at +500G to white at -500G), while panels on the right show the radial field at a height of 0.6R* above the surface (greyscale 
extends from black at +30G to white at -30G). The crosses mark the positions of the footpoints of the open field lines. 



2 EXTRAPOLATING THE CORONAL FIELD 

We write the magnetic field B in term of a flux function 
^ such that B = — V\P and the condition that the field 
is potential (V X B — 0) is satisfied automatically. The 
condition that the field is divergence-free then reduces to 
Laplace's equation V 2 *]/ = 0. A solution in terms of spheri- 
cal harmonics can then be found: 



/J [ai m r l + bin 



-O+i)' 



(i) 



l — 1 m— — l 



where the associated Legendre functions are denoted by P; m . 
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The coefficients ai m and bi m are determined by imposing 
the radial field at the surface and by assuming that at some 
height R s above th e surface the field becomes rad ial and 
hence B e (R s ) = (|Altschuler fe Newkirk, Jr. 196S|). Since 



large slingshot prominences are observed on AB Dor mainly 
around the co-rotation radius which lies at 2.7J?* from the 
rotation axis, we know that much of the corona is closed out 
to those heights and so we set the value of R a to 3.4i?*. In 



1 



order to calculate the field we used a code originally devel- 
oped by van Ballegooijen et al. (1998). 

We use three data sets, obtained on 1995 Dec 11-13, 
1996 Dec 23-25 and 1998 January 10-15 (see Fig. |l|). In each 
case we extrapolate the coronal field and distinguish between 
those field lines that are closed and those that are open. Fig. 

shows some sample open field lines and also the separator 
surfaces that separate the regions of open and closed field. 
In all three cases the global topology is similar and so we 
show only the results for the 1995 data set. There are two 
dominant regions of open field lines formed at mid-latitudes 
about 180° of longitude apart. They are of opposite polarity 
and hence form large helmet streamers, below which the field 
is predominantly closed. 



3 THE EFFECT OF THE UNOBSERVABLE 
HEMISPHERE 

Since the rotation axis of AB Dor is inclined at some 60° to 
the observer, only one hemisphere can be imaged reliably. 
The global structure of the coronal field, however, depends 
on the way in which field lines originating in the observable 
hemisphere connect to the hidden hemisphere. Close to the 
surface the effect of this missing information is negligible as 
the small-scale fieldlines connect locally to the surface. The 
large-scale field is much more likely to be affected. While 
we have no way of determining the field structure in the 
hidden hemisphere, we can assess the extent of its influence. 
We do this by creating an artificial surface map in which 
the 1995 data set forms one hemisphere and the 1996 map 
forms the other. This produces a global field structure where 
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Figure 4. Global field topology for a symmetric combination of the 1995 and 1996 surface maps where positive polarity open regions 
are at the same longitude in each hemisphere ( see a lso F ig, h h. The top three panels show a view from longitude 220° while the bottom 
panels show a vi ew fr om longitude 130°. Panels 4(a) and k(djshow the separator surf aces that sepa rate the closed and open field regions. 
Panels |4(b)| and 4(e)| show the closed field lines just inside the separator and panels 4(c) and 4(f) show the open field lines just outside 



it. In panels 4(c) 



and 4(f) the surface radial map has been painted onto the stellar surface. 



the two hemispheres have similar large-scale field topologies. 
It is of course possible that the invisible hemisphere has a 
different structure to the one that we can observe. If, for 
example, the dynamo were to excite a mixture of modes 
that are symmetric and antisymmetric about the equator 
then the magnetic structure of the two hemispheres could 
be quite different. In the light of our limited information, 
however, we choose to make the simplest assumption, that 
the lowest-order field in the two hemispheres is either purely 
symmetric or antisymmetric about the equator. We obtain 
these two cases by selecting the alignment of the maps for 
the two hemispheres. In one case we choose a symmetric 
alignment, such that the two positive regions of open field in 
each hemisphere are at the same longitude. In the other case 
we place them 180° of longitude apart (see Fig. ||). In the 
"symmetric" case, we therefore have very extended coronal 
holes that reach from one hemisphere into the other, while 
in the "anti-symmetric" case, the low-latitude field regions 
that were originally open now connect across the equator to 
form closed field regions. 

In the symmetric case (Fig. ^), the closed corona is 
confined to a torus that covers both poles and has an axis 
that connects the longitudes of the open field regions. The 
open field regions extend across both hemispheres, reach- 
ing to fairly high latitudes, but not to the poles. In the 
anti-symmetric case (Fig. ^), the closed corona is a complex 



surface that appears to be the sum of two tori: one torus is 
similar to that for the symmetric case and the other torus 
lies in the equatorial plane and has its axis parallel to the 
rotation axis. This second torus is formed by field lines from 
low to intermediate latitudes in each hemisphere connecting 
across the equator. 

In both cases, the structure of the high- latitude field and 
the longitudes of the coronal holes in the visible hemisphere 
are very similar. Indeed, they are very much the same as in 
the case where there is no field added into the hidden hemi- 
sphere at all. The main difference is in the structure of the 
low to intermediate latitude field. In the symmetric case, any 
X-ray emission would have to come from two distinct lon- 
gitude bands (indeed, any prominences formed would have 
to originate within these bands too). This is in conflict with 
the wide range of rotation phases at which prominences are 
observed. The antisymmetric case would allow prominence 
formation at any longitude and also perhaps a greater X- 
ray emission measure, since a greater volume of the corona 
is closed. In the symmetric case, the volume filling factor 
(i.e. the volume of the closed corona as a fraction of the en- 
tire volume out to the source surface) is 0.24, whereas in the 
antisymmetric case it is 0.38. 
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Figure 5. As Fig. |l| except that the alignment of the two maps is such that both hemispheres are anti-symmetric about the equator, 
with positive polarity open regions 180° apart in longitude in the two hemispheres (see also Fig. The top three panels show a view 
from longitude 308° while the bottom panels show a view from longitude 128°. 



4 THE EFFECT OF FLUX HIDDEN IN THE 
DARK POLAR CAP 

Doppler images of AB Dor consistently show that the polar 
regions above latitude « 70° — 80° are dark. The appearance 
of a dark area suggests that the field in this region is strong 
enough to inhibit convection, but Zeeman-Doppler imaging 
recovers little if any field here. This is principally because 
the Zeeman signal is suppressed in areas of the surface that 
are dark. We are therefore unable to determine either the 
polarity or strength of this field (or indeed to determine if 
it is of uniform or mixed polarity) . 

While a dark polar cap of mixed polarity would have 
limited effect on the global field, the same is not true of a 
unipolar region. In the limiting case where this polar field 
dominates the field structure, the closed corona would be 
in the form of a torus lying in the equatorial plane and 
the closed field regions seen at the pole in Figs. ^ and ^ 
would be replaced by open field. In order to determine the 
extent to which such a polar field might affect the global field 
structure, we have added dipolar fields of different strengths 
to the observed surface map for 1995. The surfaces bounding 
the closed corona are shown in Fig. []. A dipole field whose 
strength at the pole is only around 100G is sufficient to 
sweep the polar regions clear of closed field lines. Such a 



field could easily pass undetected in the Zeeman-Doppler 
maps, but would be insufficient to suppress convection. 

While the volume filling factor (i.e. the volume of the 
closed corona as a fraction of the entire volume out to the 
source surface) of all of these models is almost identical at 
0.25, what distinguishes them is the orientation of the torus 
that forms the boundary of the closed corona. In particular, 
there are significant differences in the latitudes at which the 
open field is found (and hence from which a wind could 
escape). The ratio of the open flux to the total flux is in 
fact a maximum when the imposed polar field dominates 
over the observed field. It should be stressed that although 
we have shown results for modelling the polar cap by the 
addition of a dipolar field to the observed field, the addition 
of polar spots containing the same flux as the dipole case 
gives qualitatively the same results. 



5 THE LIMITING CASE OF A POLAR FIELD 
DOMINATING THE GLOBAL TOPOLOGY 

The inevitable result of imposing a polar field that is strong 
enough to suppress convection is that this field dominates 
on the largest scales. The global field topology is just that of 
a dipolar field with a source surface imposed. In particular, 
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Figure 6. Surfa ces se parati ng cl osed an d op en fiel d reg ions for differe nt dip olar fie lds ad ded to the surface maps. Field strengths added 



are: 3(a) 1000G, 3(b) 200G, S(c 



drawn in each case. 



100G, 



6(d) 



20G, 



6(c)| 0G, |6(fJ -20G, |3(g)| -100G, |6(h) 



-200G. 



6(i) -1000G. Sample open field lines are 



as can be seen from 3 (a 



or 3(i), the polar regions are open 



and the polar hole extends down to a latitude that depends 
on our choice of the source surface position. This result ap- 
pears to be in confli ct with the BeppoSAX observations of 
Maggio et al. (2000) that suggest the presence of closed loop 



structures at high latitudes. 

In order to explore the implications of such a polar field, 
we can examine the tractable case of a purely dipolar field. 
We therefore take the I = 1, m = component of (^) and 
impose the boundary conditions 



B r (r : 
B e (r- 
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(6) 



where M is the dipole moment for a purely dipolar field 
and can be defined as M = B r (r = i?*,0 = 0)i^/2. These 
boundary conditions give a "pseudo-dipole" field which can 
be written in terms of a correction to the classical dipole 
field which allows for the effect of the source surface: 



2M cos 6 



M sin 



r 3 + 2R 3 S 
Rl + 2R 3 S 
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(7) 
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The equation of a field line is sin 2 9 = Ar/(r 3 + 2R 3 S ). The 
last closed field line passes through 8 — n /2 at r — R s and 
so has A — 3Rg. It connects to the stellar surface at a co- 
latitude Oo where 



! e = 



3R S R* 
R 3 + 2R 3 ' 



(9) 



Fig [?] shows how the latitude of this last closed field 
line varies with the imposed source surface radius R s . For 
values of R 3 below about 3R* this limiting latitude decreases 
rapidly, and the size of the polar hole increases accordingly. 
The source surface would have to be placed beyond 6R+ to 
push the polar hole above 60° latitude and to allow part of 
the cl osed field region to b e uneclipsed, as in the Beppo-SAX 
flare (Maggio et al. 2000). This would reduce the fraction 
of open flux, since the flux of open field through the stellar 
surface in the upper hemisphere is just 2nRl J®° B r sin 6d9. 
Hence the ratio of this open flux to the total flux through the 
upper hemisphere is simply sin 2 0o. As the source surface 
is moved further from the stellar surface, the fraction of the 
flux that is open decreases. 

While the amount of open flux is important to the rate 
at which angular momentum can be lost in a stellar wind, 
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Figure 7. Latitude of the boundary of the polar hole as a function 
of the source surface radius. Beyond this spherical surface the field 
lines are forced to be open. 
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Figure 8. Volume filling factor of the closed corona as a function 
of the source surface where the field lines are forced to be open. 



it is the volume of the closed corona that is relevant to the 
amount of X-ray emission that is produced. For this pseudo- 
dipole field the volume filling factor can be calculated by 
integrating the volume under the last closed field line whose 
path 6 = h(r) is defined by sin 2 9 = 3R 2 r/(r 3 + 2R 3 S ): 



Volume = 
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r 
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The dependence of this volume on the position of the source 
surface is shown in Fig. |^. As the source surface moves closer 
to the star, the filling factor drops rapidly. While the X-ray 
emission measure also depends on the density, this shrink- 
ing of the available volume will have an increasing effect as 
the source surface approaches the stellar surface. Jardine & 
Unruh (1999) have shown that if rapid rotation strips open 



polar field dominates, the only way to have the corona ex- 
tending to latitudes above about 60° is to have the source 
surface beyond about 6i?*. This is consistent with the obser- 
vations of prominences forming between 3 and 5 i?*, but does 
raise the problem of containment. The equatorial co- rotation 
radius on AB Dor is at some 2.7R* from the rotation axis. 
Beyond this point, for an isothermal corona, the density and 
pressure of the corona start to rise. The magnetic pressure 
falls with height, however, and so inevitably at some point 
the plasma pressure will exceed the magnetic pressure. The 
plasma will be able to distort and ultimately open up the 
field lines. Our imposition of a source surface is a crude way 
of modelling this process. If we place the source surface far 
beyond the co-rotation radius, an implausibly strong field is 
required to confine the plasma. 

We can quantify this problem by determining the 
plasma pressure as a function of height for the case where 
the polar field dominates. If we impose hydrostatic equilib- 
rium along a field line, then the pressure is simply 



(12) 



p = po exp ( — (p g s ds 



Here T is the temperature, g s is the component of the effec- 
tive gravity along the field line, i.e. g 3 = (g.B)/|B| and 

g(r,0) = (-GM+/T 2 +aAsin 2 

where u> is the angular velocity. Using 
as 



ui r sin9 cos9j , (13) 
we can write this 
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where $ g and $ c are the surface ratios of the gravitational 
and centrifugal energies to the thermal energies, i.e. 

GhLjR* 



4> 



kT/m 
lu 2 RI/2 
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The constant C varies from one field line to the next. If we 
focus on the equatorial plane, where each field line has its 
maximum extent r m , then C is given by 



C 



2R+(r 3 m + 2R 3 ) 



r m (Rl+2Rl)- (15) 

The variation with height of the magnetic pressure is 
more straightforward. Looking just at the equatorial plane, 
we find that 

2M (~r 3 + R 3 S ) 



\B\ = 



(m + 2Ri)- 



(16) 



If we now take the ratio f3 of the plasma pressure to the 
magnetic pressure (B /2/x) we find that it rises rapidly with 
distance above the stellar surface. If, as in Fig. ^| we put the 
source surface at 6R t , we find that the ratio reaches unity 
at around r — 3.3i?* for a base pressure of lPa, a tem- 
perature of 10 7 K and a field strength at the pole of lkG. 



This pressure is consistent with the results of Maggio et al 



the out 3r parts of stellar coronae then this shrinking of the (2000) . The height at which j3 = 1 is in fact fairly insensi- 



coronal volume with increasing rotation rate could explain 
the observed saturation and supersaturation of the X-ray 
emission . 

What is clear from Fig. Q is that in a star where the 



tive to our choice of parameters. Placing the source surface 
further out (say to 10i?*) or dropping the coronal tempera- 
ture from 10 7 K to 10 6 K moves this point out closer to 4R*, 
but in either case, the plasma pressure is one or two orders 
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by having a smaller polar spot surrounded by a ring of op- 
posite polarity as in the models of Schrijver & Title (2001). 
Field lines in the polar spot will connect preferentially to 
the surrounding opposite-polarity ring, forming closed field 
regions at high latitudes. 

If indeed young magnetically active stars do have 
mixed-polarity regions in their dark polar caps, there are 
implications for models of pre-main sequence stars. Many 
earlier models for the structure and formation of accretion 
disks and jets have rested upon the assumption that the 
star itself has a dipole-like field (|Shu et al. 1994 Lovelace 



Romanova fc Bisnovatyi-Kogan 1995 ; Miller fc Stone 1997 ) 



12 3 4 

Radius (R.) 

Figure 9. The ratio /3 of plasma to magnetic pressure along a 
field line as a function of radius. The field lines are forced to be 
open at a radius R s = 6i?* . 

of mag nitude greater than the magnetic pressure before the 
source surface is reached. Even with a base pressure as low 



as O.OlPa with the source surface at 10i?*, the field lines on 
which j3 > 1 do not extend above latitude 60°. 



6 CONCLUSIONS 

We have used Zeeman-Doppler images of the surface mag- 
netic field of AB Dor to extrapolate the coronal magnetic 
field. This field has a clear non-axisymetric structure that 
becomes very apparent upon tracing the open field lines 
that fill much of the coronal volume. These open field lines 
originate in two opposite-polarity mid-latitude regions, sepa- 
rated by about 180° of longitude. The surface separating the 
closed and open field regions forms a torus that passes over 
the visible pole. We have tried to allow for the two regions 
in which flux could be missing from these maps: both in the 
unobservable hemisphere and in the dark polar cap. Allow- 
ing for field in the unobservable hemisphere changes the field 
line connections at low latitudes, but the high-latitude field 
topology is largely unaffected. If we allow for the flux that 
may be concealed in the dark polar cap, we find that the po- 
lar field lines become open with the addition of a polar field 
so weak that it would not suppress convection sufficiently 
to cause the polar regions to be dark. If we impose a polar 
field strong enough to give a dark polar cap, we find that the 
large scale field is similar to that of a dipole with a source 
surface imposed. 

This "pseudo-dipole" field does not however explain 
the BeppoSAX flare observations which imply high latitude 
closed loops. The boundary of the polar hole reaches down 
to latitudes well below that where observations suggest the 
footpoints of the flaring loop must have been located. The 
only way to reduce the size of the polar hole sufficiently is to 
move the source surface out to beyond 6i?*. This, however, 
requires the pressure in the tops of the largest loops to be 
much greater than the magnetic pressure. The only reason- 
able solution seems to be to allow the polar flux to be of 
mixed polarity. This could be achieved by having the dark 
polar cap composed of many smaller spots of mixed polarity, 
as in the flux-emergence models of Schiissler et al. (1996), or 



More recently, Agapitou & Papaloizou (2000) and Terquem 
& Papaloizou (2000) have shown that departures from a po- 
tential field or a mis-alignment of the dipole field with the 
stellar rotation axis (i.e. a non-axisymmetric structure) can 
influence the angular momentum transfer between the star 
and the disk or cause an observable warping of the disk. 
Angular momentum loss in a stellar wind is also sensitively 
dependent on the surface positions from which open field 
lines escape (Solanki, S.K., Motamen, S. & Keppens, R 
). In Fig 



1997 



the number of open field lines drawn in each 
panel is directly proportional to the surface area of open field 
lines. The nature of the field in the polar cap clearly has a 
significant influence on the structure of the open field and 
hence potentially on the spin-down rate of young stars. Our 
Zeeman-Doppler images imply that allowing for departures 
from dipolar fields is certainly justified. 
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